© 2003 by the authors 

Adv. Theor. Math. Phys. 7 (2003) 667-710 



Existence of Atoms and Molecules 
in Non-Relativistic Quantum 
Electrodynamics 



Elliott H. Lieb 1 and Michael Loss 2 



1 Departments of Physics and Mathematics, Jadwin Hall 
Princeton University, P. O. Box 708, Princeton, NJ 08544 

lieb@princeton.edu 

2 School of Mathematics, Georgia Tech, Atlanta, GA 30332 
loss@math.gatech.edu 

Abstract 

We show that the Hamiltonian describing N nonrelativistic elec- 
trons with spin, interacting with the quantized radiation field and 
several fixed nuclei with total charge Z, has a ground state when 
N < Z+l. The result holds for any value of the fine structure constant 
a and for any value of the ultraviolet cutoff A on the radiation field. 
There is no infrared cutoff. The basic mathematical ingredient in our 
proof is a novel localization of the electromagnetic field in such a way 
that the errors in the energy are of smaller order than 1/ L, where L is 
the localization radius. 



e-print archive: http://lanl.arXiv.org/abs/math-ph/0307046 

1 Work partially supported by U.S. National Science Foundation grant PHY 01-39984. 
2 Work partially supported by U.S. National Science Foundation grant DMS 03-00349. 
This paper may be reproduced, in its entirety, for non-commercial purposes. 



668 



EXISTENCE OF ATOMS IN NON-RELATIVISTIC QED 



1 Introduction 

The existence of atoms and molecules in the framework of the Schrodinger 
equation was proved by Zhislin [14] for fixed nuclei when N < Z + 1. That 
is to say, the bottom of the spectrum of the TV-electron Hamiltonian is a 
genuine iV-particle bound state that satisfies Schrodinger's equation with 
some energy E, for each choice of the locations of the nuclei. (Here N is the 
number of electrons, each of charge — e, and Ze is the total charge of one 
or more fixed, positively charged nuclei.) The main physical result of the 
present paper is the proof of the same thing when account is taken of the 
ever-present quantized electromagnetic field. The interaction of this field 
with the electrons (but not the field itself) necessarily has an ultraviolet 
cutoff \k\ < A (in order to have finite quantities), but we emphasize that no 
infrared cutoff is used here. 

If the fine structure constant a = e 2 /hc and A are small enough, the 
result follows from [1], but our result holds for all values of these parameters. 
In a recent paper Barbaroux, Chen and Vugalter [3] developed a new method 
that shows the existence of ground states for two-electron molecules with 

2 < Z + 1 (e.g., the Helium atom). Although they do not have to require 
that the perturbation is small when compared to the ionization energy as in 
[1], they have to impose restrictions on the various parameters since their 
works relies on the existence of the zero momentum ground state of the 
Hamiltonian of an electron interacting only with the radiation field. This 
has been established in [4] but only for sufficiently small coupling constants. 
The method of [3] is different from ours. 

Our work (and [3]) relies on earlier work with Griesemer [7] where it was 
shown that a ground state exists provided a "binding condition" is satisfied, 
and it is this condition that is proved in [3] for the restricted N = 2 case 
and for the general case here for TV < Z + 1. If E v (N) denotes the bottom 
of the spectrum of the Hamiltonian H V (N), which includes the Coulomb 
attraction of the electrons to the fixed nuclei of various positive charges 
Z\e, . . . , Zk& with Z = Yl Zj, an d if E°(N) denotes the bottom of the 
spectrum of H°(N) — the "free-electron" Hamiltonian in which there are 
no nuclei, but the electron-electron Coulomb repulsion is included — then 
the binding condition is 

E v (N) < min {E v (N') + E°(N - N') : < N' < N} . (1.1) 

This binding condition, incidentally, is the same condition that Zhislin de- 
rived for the Schrodinger equation without the quantized electromagnetic 
field, and which he verified for N < Z + 1. 
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The inclusion of the quantized electromagnetic field presents two main 
difficulties. One is that if the bottom of the spectrum contains an eigenvalue 
it is not an isolated eigenvalue, as it was in [14]. Rather, the bottom of the 
spectrum is always the bottom of the essential spectrum because one can 
create arbitrarily many, arbitrarily soft photons. It is not easy to find an 
eigenvalue when it lies in the continuum. This problem was solved in [7] 
under condition (1.1). 

The second main problem, which complicates the proof of (1.1), comes 
from the fact that each electron carries a virtual cloud of photons. This 
cloud may have substantial energy and when two electrons are near each 
other (whether bound or not) the interference of the photon clouds must be 
taken into account. In general, this is a highly non-perturbative effect. Our 
way around this difficulty is to prove that the photon clouds can be localized 
(i.e., effectively eliminated outside a ball of radius L surrounding the electron 
or the atom) in such a way that the error induced in the energy of the cloud 
is smaller than L~( 1+£ \ and thus the direct Coulomb interaction, which goes 
as L -1 , is dominant — as it was in the original paper [14]. A closely related 
effect is that even in the absence of an external potential electrons interact 
with each other. In such a case their dynamics is governed by H°, which 
contains the electron-electron Coulomb repulsion. Nevertheless, it is not 
inconceivable that the quantized field, which interacts simultaneously with 
all the electrons, might cause binding among the "free" electrons. While this 
is unlikely it has never been disproved and we must not assume in our proof 
that E°(N) = NE°(1). 

We are grateful to the anonymous referee who made many valuable sug- 
gestions and who helped us understand some conceptual matters about pho- 
tons. 



2 Basic Definitions and Concepts 

The Hamiltonian under consideration, in appropriate units, is the Pauli-Fierz 
Hamiltonian and is given by 



A? 



H V (N) = [iPj + V^A( Xj )f + |yW, • B( Xj ) + V(xj) 



3=1 




(2.1) 
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Here, g is some constant (close to 2, physically) and the vector aj is the set 
of three Pauli spin matrices for electron j. (Owing to the ultraviolet cutoff 
there is no restriction to \g\ < 2, as there would be without a cutoff [12].) 
The operator pj denotes — iV acting on the coordinate of the j-th electron. 

The potential V is the potential of K > 1 nuclei with positive charges 
Zi,..., Zk and locations Ri, . . . , Rk € K 3 - 

K 

V(x) = -J2Zj\x-R j \- 1 . (2.2) 



Remark: The truth of our main theorem (3.1) — and its proof — does 
not require that V{x) be given by (2.2). In addition to the general condition 
[7, eq. (5)], we need only the condition that there is some radius p such 
that (V(x)) < —Z/\x\ for \x\ > p, where ( • ) denotes spherical average. 
Similarly, the repulsion | X% Xj I Cclll be replaced by W(xi — Xj) provided 
(W(x)) < l/\x\ for all \x\ > p. 

The free Hamiltonian H°(N) is similar to H V (N), but without the at- 
traction to the nuclei, i.e., 

N 

H°(N) = Y^ [(Pi + VSA(xO) 2 + ^<Ti-B{xij\ | ._ i +H f . 

1=1 i<3 J 

(2-3) 

Note that the Coulomb repulsion among the electrons is included. The 
reason for including the Coulomb repulsion is (as stated above) that we do 
not know whether the electrons bind to each other through the interaction 
with the electromagnetic field, i.e., the electrons may not separate in the 
lowest energy state (if there is one) . 

The (ultraviolet cutoff) magnetic vector potential is defined by 

A W = 7 tPj 7W 2a W + «a(%-* X ) dk , (2.4) 

where the function Xa 1S a smooth, radial function in k space, that vanishes 
outside the ball whose radius is the ultraviolet cutoff A. We denoted the cre- 
ation and destruction operators of photons of momentum k and polarization 
A by a\(k) and a*Ak). This unusual notation is used since we shall later in- 
troduce the creation and destruction operators in configuration space, a\(y) 
and a* x (y), which act on the Fourier transformed functions in Fock space. 
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The magnetic field is B(x) = curL4(x). The operators a x ,a* x satisfy the 
usual commutation relations 

[ax(k),at(q)]=8(k-q)8x,u , [a x (k),a v (q)] = 0, etc (2.5) 

and the vectors S\(k) are the two possible orthonormal polarization vectors 
perpendicular to k and to each other. 

The vectors £\(k) have to be discontinuous functions of k on every sphere 
of fixed | A; | -value because it is not possible to "comb the hair on a sphere". 
However, the only physical quantity, 

E4(*)4(*) = *j-^> (2-6) 

A=l 1 1 

is discontinuous only at the point k = 0. For the rest of this paper we choose 
the polarizations vectors to be 



ei(*) 



(fc 2 ,-fci,o) 



e 2 (k) = j£|Aei(fc). (2.7) 



Let us emphasize here that some smoothness of the function x A is essen- 
tial for our arguments since this guarantees that the coupling functions 

h\(y) = i- J ^M £ \(k)e- ik -ydk (2.8) 



has a suitable decay as \y\ — ► oo. If we did not have the discontinuous 
function e\(k) in (2.8) then h(y) would decay as |y| _5//2 as \y\ — > oo. (Proof: 
|fc| -1 / 2 is the Fourier transform of |y|~ 5 / 2 in the sense of distributions [10, 
Theorem 5.9]. The Fourier transform of xa is re& l analytic and decays faster 
than any inverse power of \y\. Hence, the convolution of x with |y|~ 5 / 2 
decays like |y| -5 ^ 2 . With a sharp cutoff it would decay only like \y\~ 2 which 
turns out to be insufficient for a good localization of the photon states. 

This analysis of h shows that we have to be circumspect about the choice 
of the polarization vectors. Their discontinuity will spoil the |y|~ 5 / 2 decay, 
but it is important to get better decay than |y| -2 . In Lemma B.l of Appendix 
B it is shown that with our choice (2.7) of the polarization vectors the 
coupling functions have sufficient decay in the sense that f \y\ 2l \h\{y)\ 2 dy is 
finite for all 7 < 1. Thus, in an average sense, the coupling functions decay 
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almost as fast as \y\ 5 / 2 . We made no attempt to optimize the choice of the 
polarization vectors. 

While polarization is physically measurable, the polarization vectors are 
not. They are merely a basis. It is odd, therefore, that their mathematical 
definition plays a role in the spatial localization of the photon field that we 
shall construct, and which is central to our proof of the binding condition. 
It would be better to start with a formalism that contains only "divergence- 
free" vector fields as the dynamical variables instead of trying to define them 
with the aid of unphysical polarization vectors. In particular, the Fock space 
would be built over the L 2 -space of divergence-free vector fields instead of 
L 2 (g)C 2 . We shall not explore this here, but we mention that the localization 
of a divergence-free vector field, which preserves the divergence-free property 
is also a subtle matter. 

The field energy, Hf, sometimes called dT(uj), is given by 

H f=Y\ I \k\a*\(k)ax(k)dk (2.9) 

A=l,2 3 

There is no cutoff in Hf. The energy of a photon is \k\. 

Another unbounded operator of interest is the number operator 

M = V / a\{k)a x {k)dk . (2.10) 

The physical Hilbert space for this system is given by 

H(N) = A N L 2 (R 3 ;C 2 ) ® JF (2.11) 

where the wedge indicates that the electron wave functions are antisymmetric 
under the exchange of the particle labels. Thus, the functions in the space 
TL(N) obey the Pauli exclusion principle. The photon Fock space is T . We 
denote the inner product of two states * and $ in the space 7i{N) or in 
Fock space alone by 

(*,$) and (*,$) , (2.12) 

respectively. If \l/ and $ are in H(N) then (\1/,<J>) makes sense and defines 
a summable function of x\, s\, ... ,xn,sn, where the space-spin 

variables of the j-th electron. 

It is desirable that the above Hamiltonians be selfadjoint on certain do- 
mains and this has been worked out, e.g., in [9]. In this paper we will always 
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be talking about the Friedrichs extension of the symmetric operators H#(N) 
(where # is or V). The form domain will consist of all states for which 
each term in the operators has a finite expectation value. Accordingly, we 
define the ground state energy E# (N) for the Hamiltonian H# (N) by 

E*(N) = inf | H*(N)^ : * G H(N) , \\V\\ = l} . (2.13) 

The numbers E#(N) are finite. This follows from Lemma A. 4 in [7] together 
with the fact that the Coulomb potential is form bounded with respect to 
p 2 = -A. 

A few remarks concerning the Fock space T are in order. It is built over 
the space L 2 (M 3 )<g>C 2 ; the second factor takes into account the polarizations. 
Let {fi}, i = 1, 2, . . . , be an orthonormal basis for L 2 (M 3 )<g>C 2 . Then, vectors 
of the form 

\h,m i; ... ;i n ,m n ) = ^=^=a*(/ n r i ---a*(/ l J m "|0) , (2.14) 

\fm\\ ■ ■ ■ m n \ 

constitute an orthonormal basis for J 7 , the occupation number basis. In 
(2.14) n is an arbitrary nonnegative integer (with n = denoting the vac- 
uum vector |0)) , the indices ii, ■ ■ ■ ,i n are all different, the m-i are all positive 
integers, a*(f) is an abbreviation for J2\ a *\(f>>) an d fx = f{k, A) is a func- 
tion in L 2 . Thus, any state $ G T can be uniquely written as 

,mi; ... ;i„,m n \h,mi] ■■■ ',in,m n ) , (2.15) 

where the n = term in (2.15) is just </>o|0) with 4>q G C. The inner product 
is given by 

= E E E I0u,mi; ... ;wnJ 2 ■ (2-16) 

n>0 ii<i2< ■■■ <i n ■■■ ,m in 

This representation has the advantage that the symmetry in the photon 
variables is automatically taken care of. It is particularly useful when dealing 
with product states. Consider a state <5 whose photons are all localized in 
a closed region y C M 3 . This means that all the fi(y, A) appearing in (2.14) 
and in (2.15) vanish if y £ y. Likewise, consider a state ^ whose photons 
are all localized in a closed region Z C M 3 which is disjoint from y. Pick 
an orthonormal basis {fk} in L 2 {y) (g> C 2 and an orthonormal basis {g{\ in 
L 2 (i?)®C 2 . Clearly, the two algebras of creation and annihilation operators 
generated by a*(/fc) and a^{g^) commute. If 

$ = E E E <f>h,PU - ;in, P n \h,Pi; ... ;i n ,Pn)y (2.17) 

n>0 h<l2< ■■■ <in Pi, ••• ,Pn 
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and 

* = E Yl Yl ^hm; - ;j k ,q k \h,Qi\ ■■■ \3k,Q.k)z ( 2 - 18 ) 

then we define the product state E by 



)y® \h,qi; ■ ■ ;jk,gk)z 

(2.19) 

where 

In, pi; ••• ;im,Pm)y®\ji,qi; ■■■ ;jk,qk)z = 



(2.20) 



By a simple calculation we find that 

(S,S) = . (2.21) 
Further, if / is a function supported in y, then 

(E, o*(/)a(/)H) = <<&, a*(/)a(/)d>) *> . (2.22) 
Likewise, if / is supported in y and 5 in Z, then 

(~,a*(f)a(g)E) = ($,«*(/)$) <\MG/)*> . (2.23) 

Quite generally, we have, for normal-ordered, bilinear expressions, the 
following formulas (in which (3, 7 denote linear forms in the annihilation 
operators a, and hence (3*, 7* are linear forms in the creation operators): 



<H,/3 7 H) = (tf,/37*> + 7^) 

+ /3 *) ($, 7 $) + 7 *) ($, /3 $) 
(S, /?* 7* 5) = (*, (3* 7* *) $) + *) <<*>, /?* 7* $ ) 

+ (*, /?* *) <$, 7* $) + (*, 7* *) <$, /?* $) 
(S, /F 7 H) = (*, /T 7 *) + *> <<M* 7 *) 

+ (*, (3* *) ($, 7 $) + (*, 7 *) ($, /?* $) (2.24) 
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A formula of the type (2.24) does not exist for anti-normal-ordered prod- 
ucts /?7*. We shall have no need of such terms, however, because the only- 
source of such terms is A{xi) 2 in the electron kinetic energy. If we denote 
the part of (2.4) coming from a\(k) by (3(x) and the remainder by (3*(x) (see 
(2.29)) then 

A(x) 2 = P{xf + (3*{x) 2 + 2p*(x)p(x) + C (2.25) 

where 

Thus, apart from a fixed, finite number aNC, which is strictly proportional 
to N, (and which is, therefore, independent of any decomposition of the 
system into clusters) we can (and henceforth shall) replace A(xi) 2 by the 
normal-ordered 



: A{ Xi ) 2 := f3{ Xi ) 2 + (3*{ Xi ) 2 + 2(3* (xi)f3(xi) . (2.27) 

Formulas (2.24) continue to hold for vectors $ in the physical Hilbert 
space 7i(N), with the replacement of ( , ) by ( , ). In this case, the coeffi- 
cients, <f>i 1>mi; ... ;i n ,m n , are (antisymmetric) functions of the electron space- 
spin coordinates, Xi, Sj. 

It is convenient to introduce the operators given by 

ax(y) = 7^3 / Mk)e tk - y dk . (2.28) 
Then the vector potential can be written as 

2 

A\x) = J2ax(h\(x - •)) + a\{h\{x - ■)) . (2.29) 

A=l 

The action of the operators a\(K l x (x — •)) is given by 

[a\(h\(x- -))^]n(yi,Ai;--- ;y n ,K) 

= + 1 J h\(x-y)[9] n+1 (y,X; yi ,Xi;--- ;y n ,K)dy ■ (2.30) 

A convenient expression for a\(h\(x — ■)) is the formula 

a x (h\(x -0) = / a x (y)h\(x-y)dy . (2.31) 
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The number operator and the field energy can be expressed in terms of 
the operators a\(x) by 

2 r 

(d>,AA<I>) = (2tt) 3 £ / \\a x (x)*fdx , (2.32) 

A=l J 

and 

2 

($, H f S) = (2tt) 3 ( o a0)*, \/— A a A (-)$) (2.33) 

A=l 

which by eq. 7.12(4) in [10] can be rewritten as 

47r z^y dxd v ■ ( 2 - 34 ) 

By the previous considerations we have for the product state S that 
E(a A (x)S,a A (y)E) = 

A 

A i 
+ J](a A (x)*,*)(*,a A ( J /)*) + j;(*,a A (y)*)(a A (a:)*,*) , (2.35) 

A A 

and hence we obtain for the field energy of E the expression 



(E, H f Z) = (<D, H f $) (tf , *) + $) (*, if,*) 

(a A (x)$, $) (*, a A (y)tf ) + ($, a A (y)$) (a A (x)*, *) 



- 8vr £ K y 



|x — y| 4 



(2.36) 



The x integration in the first term of the last integral runs over the set 
y while the y integration runs over the set Z and similarly in the second 
term the x integration runs over the set Z while the y integration runs over 
the set y. Hence the last expression is well defined as long as the distance 
between the sets y and Z is positive. This term expresses the fact that the 
field energy is a nonlocal operator and this nonlocality is one of the main 
obstacles to be overcome. 

In general, the states $ and will depend on the position and spin 
variables of the various electrons and hence the product state (2.19) has to 
be antisymmetrized over the electron labels. It is straightforward to check 
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that the expression (2.36) continues to hold also for such states. (When 
different groups of electrons are involved an antisymmetrization is required, 
however, as discussed in (3.8).) 

We need one more concept before stating our main theorem. It will 
be necessary to localize both the electrons and the photon field. As far 
as the electrons are concerned it is useful to define what we mean by a 
symmetrized product of n domains in M 3 . If B±, . . . B n are n domains 
(open sets) in M 3 then the symmetrized product, Q, is a domain in (IR 3 )™ 
given by 

n(Bi, ...,B n )= |J B vl x B„ 2 x ■ ■ ■ x B m , (2.37) 

where S n is the group of permutations of n labels. It might be useful to 
illustrate this when n = 2. Then we have B x C IR 3 , B 2 C M 3 and Q(Bi,B 2 ) = 
(B l x B 2 )U{B 2 x B{). This is different from (BiUB 2 ) x (B 2 UBi). Physically 
it means there is one particle in each domain f?j, but the label of the particle 
is indeterminate. If the domains overlap there may be several particles in 
one domain, of course. 

3 The Main Theorem 

The following is our main theorem. The proof given in this section uses 
several inequalities derived later on in this paper, but we present the proof 
now in order to make the main ideas clear without too many technicalities. 

THEOREM 3.1 (Binding in Atoms). The strict inequality (1.1) holds 
for all N < Z + 1, all g, all a and all A. In particular this implies that 
there exists a normalized ground state $(A) in 7i(N) for the Hamiltonian 
H V (N), i.e., ($(N),H(N)$(N)) = E V (N), and it satisfies H v (N)$(N) = 
E V (N)$(N). 

See the remark after eq. (2.2). 

PROOF: Our proof has three main parts. The first is the construction 
of a good trial function for N — N' (with < N' < N) localized, 'free' 
electrons and localized photons accompanying these localized electrons. The 
second part is the construction of a good trial function for N' localized 
electrons 'bound' to the given, fixed nuclei, together with localized photons. 
The third part consists in the construction of a trial function which is a 
product of these two functions and then showing that the energy is lowered 
(by a greater amount than the localization errors) because of a negative 
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Coulomb energy between the 'bound' system (consisting of electrons and 
nuclei) and the localized 'free' electrons. One difficulty in part 3 is that 
although the photons in the two regions are localized in separate regions, 
there is still a residual interaction between the two fields, given by the last 
term in (2.36), which has to be considered. This interaction comes from the 
fact that multiplication by \k\ in Fourier space is a nonlocal operation in 
position space. 

The general argument proceeds by induction. We know from [7] that one 
electron binds. Assuming that the binding condition holds for M electrons, 
all 1 < M < N — 1, we have to show that it holds for N electrons, i.e., 
E V (N) < mm{E v (N f ) + E°(N - N') : < N' < N}. Using [7, Theorem 
2.1] we may assume that the Hamiltonian H v (M) has a ground state for 
all 1 < M < N — 1. By the second part of [7, Theorem 3.1] we know 
that E V (N) < E°(N) for all N, since Z > and the attractive Coulomb 
potential is strictly negative. 



Parti. | From now on we set 

n = N — N' . 

Given < N' < N we shall construct a normalized state <fr(n) for the 
free electron Hamiltonian H°(n) with the property that the n electrons are 
localized in a symmetrized product 0, of some balls of radius Rq while the 
field is localized in balls with the same center but with radius L > 2Rq . The 
construction of <E>(n) is done in Theorem 4.3. It lies in the physical Hilbert 
space 7i(n) and has an energy given by 

($(n),ff°(n)$(n)) Cn f Rp\ , AI? m 27r 2 n 2 

(*(n),*(n)) " {») ( 1 + " 1 °^o)l) + ^ r , 

(3.1) 

for any 7 < 1, where C is some constant independent of L and Rq. (It does 
depend on 7 and on n, but n is bounded by N). 

In (3.1) the term (Cn/(L - 2i? ) 7 ) (Ro/L^) (1 + | log(AJ2o)|) comes from 
the energy needed to localize the field in n balls of radius L. The last term 
comes from the kinetic energy needed to localize n electrons in the n balls 
of radius Rq (Lemma 4.1). 



Parti. I According to the induction assumption at the beginning of this 
proof we may assume that 1 < N' < N— 1 and that the Hamiltonian H v (N r ) 
of the bound electrons has a normalized ground state F(N'). By [7, Lemma 
6.2] we know that this ground state is exponentially localized in the electron 
variables, i.e., if we denote by \X\ the quantity X^«=i \ x i\ then 

|| e 0*lr(JV')|| 2 < Cp (3.2) 
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for any f3 2 < mm{E v (N' - m) + E°(m) : < m < N'} - E V (N'). 

(Note: An error in the proof of this exponential localization in [7, Lemma 
6.2] was discovered by J-M. Barbaroux and the necessary correction was 
published in [6]. We are grateful to Prof. Barbaroux for pointing out this 
mistake to us.) 

Although it is not necessary to do so, we (strictly) localize T(N') so 
that all the electrons are in a common ball of radius Rq. Following that, 
we localize the photon field in a larger ball of radius L > Rq. The field 
localization is essential. The electron localization is not since it would be 
possible to use only the exponential decay of T(N'). The localization is done 
as follows. 

Let x < 1 be a smooth cutoff function with support in the unit ball 
centered at the origin and x( x ) = 1 f° r M < 1/2- Let 6 = \\^=\X{ x i/Ro) 
and T(N') = QF(N'). Since T(N') is a ground state, and hence satisfies the 
Schrodinger equation, we can deduce that the increase in energy due to the 
cutoff is bounded as follows. 



(r(N > ),H v (N')r(N'^j < E V {N') (T{N'),T{N'j) + N > CpeXP ^ 2 . 

(3.3) 



Inequality (3.3) follows from (3.2) and integration by parts as follows. 
With ( , ) denoting inner product in Fock space and dX denoting integration 
over the space-spin variables, we have 

f (er(N%Y,^j + iA( Xj )fer(N')) dx = 

J 3 

J ^eA^ej (T(N'),T(N')) dX 

+ 2 / S ( 0v ^ e ) • (r^O. ( v i + ^(^))r(AT')) dx 

3 

+ J e 2 (r(N'),^3 + ^3)) 2 nN'))dx . (3.4) 

3 

Since (V + iA(x)) 2 is a symmetric operator, the left side of (3.4) is real and, 
therefore, the right side must be real, too. The first term on the right side 
is real. The third term is also real because T(N') satisfies the Schrodinger 
equation, and hence (r(i\r'),Ej(Vj+iA(x j )) 2 r(iV , )> = (T(N'), (-E v (N') + 
real potentials)r(A'" / )), which is real. The middle term must, therefore, be 
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real too (when summed over all particles), and we can replace the inte- 
grand by its real part. This means that we can replace 2{T(N'), (V x . + 
iA(xj)),T(N')) by 2^(T(N'),V X T(N')} = V Xj (r(N'),T(N>)} since 
(T(N'),iA(xj)T(N')) is imaginary (because A(xj) is symmetric). 

Now, integrating by parts we can combine the second term with the 
first to yield - f(V Xj Q) 2 {T(N' ), F(N'))dxj. This is the error term (the last 

term leads to the principal term E V (N') (r(N'),F(N')j). This error term 

can be bounded by replacing Vx(%j/ Ro) by C ' /Rq times the characteristic 
function of the annulus between Rq/2 and Rq, for some constant C. But this 
characteristic function is bounded by exp(— ft Rq/2) exp(+f3\x\). Inequality 
(3.3) then follows from the exponential decay (3.2). 

Next one has to show that the error term in (3.3) is small when compared 
with ||r(iV)|| 2 . It follows from the exponential decay that 

||f(Ar')|| 2 > l-N'Cpe-P 1 * . (3.5) 

To see this note that x( x i/Ro) 2 > 1 — 9i, where g\ = 1 if \x{\ > Rq/2 
and <7j = otherwise. Then 6 > n«(l ~ 9i) ^ 1 — J2i9i- But g± < 
exp{-/3i? }exp{2/3|xi|} < exp{-(3R } exp{2/3|X|}. Therefore, 

nf(ivoii 2 > (t(n% (i-j29inN'?j 

> I- N'exp{-(3R }(T(N'), exp{-2/3|X|}r(iV / )) . 



Together, (3.3) and (3.5) imply (for exp{/3i? } > N' Cp) 

f{N>),HV{N>)f{N>)) N , c 

1 <E V {N') + 



(r(N'), r(JV')) R 2 exp{pR } - N' Cp 



2 N' 

<E v (N') + - T C l3 e~^ , (3.6) 
R 



where the last inequality holds provided that Rq is chosen such that (3Rq > 
\og{2N'Cp). 

The next step is to localize the photons in the state T(N') in a ball 
centered at the origin of radius L > Rq. This leads to a new state VP (AT') 
with 
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{p{N'),H v {N')V{N')) 
(^(N'),^(N')) 



F V (N'\ N ' Cfi CN ' R ° 

~ { ) + M exp{/3i? } - N> C p + (L - 2R )f Hi ' ( ] 

for all 7 < 1 and for Rq and L — 2Rq large enough. 

The construction of this function ^ is precisely the same as the photon 
field localization that led to the state $(n). It is carried out in Theorem 4.3. 

Thus, ty(N') is a state in which all the electrons are localized in a ball of 
radius i?o and the photons are all localized in a ball of radius L. Moreover 
the localization errors are small and given in (3.7). 



Part?,. | Now we put the pieces from Part 1 and Part 2 together and 
construct a trial function £ whose energy will be strictly below E V (N') + 
E°(n). 

As mentioned above, since H°(n) is translation invariant we can, by 
shifting, make sure that the photons in the state ^(iV') and the photons in 
the shifted state <&(n), live in disjoint sets. This will be the case when the 
smallest distance of the centers of the balls Bi,...B n with the center of the 
ball in which *B(N') lives is greater than 2L. 

Now we can form the product state E as indicated in (2.19). The 
state is symmetric in the photon variables by construction. It has to be 
antisymmetrized in the electron labels though, i.e., replace the products 

<Ph, P i;-;im, Pm iPji,qi;-;jk,<ik in ( 2 - 19 ) h Y 

c(N,N') ^U^-SWm^Cl),- ■ ■ ^r(JV-JV')) 

x 1 l } h,qi;-;jk,qk( z n(N-N'+i), - ■ ■ , z n{N)) ■ (3-8) 
where ir runs through all the permutations of N elements, 

c(N, N') = 1 (3.9) 

y/N\{N - N')\N'l 

is the normalization and Zj = (xj,Sj), the position and spin of the j-th 
electron. The expression in (3.9) is calculated by noting first that (p and tp 
are each antisymmetric in their electron coordinates and, second, that there 
is no overlap between <f> and tp because the x variables in the two functions 
have disjoint support. Informally speaking, the antisymmetrization in (3.8) 
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has no effect and could be dispensed with for all practical purposes since the 
operators V and the potentials that we consider here are local operators. If 
we the Hamiltonian contained nonlocal operators, such as the 'relativistic' 
v 7 — A then the antisymmetrization (3.8) would have a more profound effect 
- although (3.9) is still correct. 

Next, we calculate (with : : denoting normal ordering) 
(E,H(N) S) 

= ^S, ^2 [■ (Pi + V^Mxi)? ■ +| <Ti ■ B{ Xi ) + V(Xij\ 

+ ( H '5^M + ™ (310> 

in terms of the normalized 3>(n) and fy(N'). The field energy term has been 
explained previously in equation (2.36) and yields 



(E,H f E) = (*(n),H f *(n)) + (*(N'),H f *(N')) 

(a x (x)cD(n), $(n)) (V(N% ax(y)*(N>)) ^ 
\x — y| 4 



i J \x-y\ 

The Coulomb repulsion term is easily calculated to consist of three terms: 
E |^4^K n )) (3-12) 

\ N'<i<j<N \ X% X ^ J 

f \\*(N>)\\\x 1 ,..,x N ,)\\Hn)\\i(x N , +l ,..,x N ) dNx (3 M) 



i=l,j=N'+l 



Here the norm signs indicate that the norm has been taken in Fock space 
and in the spin space. 

The electron kinetic energy involves the calculation of terms of the form 
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= {*(N'),a(f) 2 *(N>)) (*(n), <&(«)) + «(/) 2 * W) (*(iV), *(iV')) 

+ (tf (W), a(/)tf (AT')) (d>(n), a(/)$(n)) . (3.15) 

Thus, we have that 

AT' 

^ (*(iV'), [: (ft + v^(x,)) 2 : +f ■ B^)] *(^')) *W) 

i=i 

+ ]T (<&(n), [: ^ + v^A(xi)) 2 : +|v^ • 5(x;)] $(n)) *(iV')) 

i=JV'+l 
AT' 

+ a Yl J {®(n),:A( Xi ) 2 :*(n)) ||^(iV')|| 2 (xi, . . . , XN')dx\ ■ ■ ■ eZxjv' 



N 

+ «E (*(iY'),:^) 2 :*(iV))l^( 
i=JV'+l ^ 



AT' 



(3.16) 

n)|| 2 (xAr/ + i, ..,XAr)dXAr/ + i • -dx^ 

(3.17) 



+ 2a Y.f mri),A(xi)<f>(n)) • (*(JV'),Pi*(JV')) (xi, • • • , xjvO^i • • • dx N , 

(3.18) 

+ 2a / (^O.^M^'))- 

i=JV'+l ^ 

($(n),pj$(ra)) (xat/ + i, . . . ,x N )dx N , +1 ■ ■ ■ dx N (3.19) 

N 

+ 2a^2 {V( N '), A(xi)V(N')) ($(n), A(x;)$(n)) dxi • • • dx N (3.20) 

AT 

+ fv^ E / (*(AO,£(x,)*(iV'))- 

i=JV'+l 



AT' 



($(n),<7i$(n)) (xjv+i, • • ■ ,XAr)dxAr/ + i • • • dx N (3.21) 



+ |v^E/ (*(n),-B(xO*(n))- 

o** (W)) (xi, . . . ,sjv/)dxi • • • dx N , 



i=i 



(3.22) 
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Finally, the last and most important term 

(~,f>0r;)~) = ^N'),Y,V( Xi MN'^ (*(n),*(n)) 

+ U(n), Y, V(xi)*(n)\ (*(N'),*{N')) . (3.23) 

V i=N'+l J 

Lemma 5.6 allows us to show that the terms (3.16, 3.17) are of order 
L~ 27 for any 7 < 1. This follows from the fact that in our trial function the 
electrons are localized in balls of radius Rq , and so the distance D in Lemma 
5.6 between any electron and the localized photon field of the subsystem to 
which the electron does not belong is at least L — Rq. Since we can easily 
choose L large and Rq to be an arbitrarily chosen small constant times L 
we conclude, from Lemma 5.6, that these terms are of order L~ 27 for any 
7 < 1. 

The terms (3.14, 3.23) taken together are the terms that will give us 
binding. We shall show that, after averaging over rotations, the two terms 
add to — {Z — N')/SL, which is less than — pos. const. /L according to our 
hypothesis. 

'Averaging over rotations' means the following. We fix the state ^(N') 
of the electrons bound to the nuclei and their field, which extends out a 
distance L from the origin. On the other hand, the state of the n unbound 
electrons was called 3>(n), but actually there are infinitely many states we 
could use. That is, we start with one $(n) and consider all rotations of 
it about the origin. The average Coulomb interaction (i.e., the average of 
(3.14) and (3.23)) is the same as if the bound electron state ^f(N'), including 
the nuclei, was rotated about the origin. However, the average potential 
generated by the latter average over rotations at a point x would be exactly 
(Z — N')/\x\ provided \x\ > L. This is Newton's theorem [10, Theorem 9.7]. 
Therefore, there exists a rotation so that the Coulomb interactions (3.14) 
and (3.23) are as if the inner state were a point charge located at the origin 
and of strength Z — N' . 

We now choose 3>(n) so that one of the balls of radius L in which the n 
electrons and the field reside is tangent to the ball of radius L in which the 
bound electrons and field reside. By averaging over rotations we may assume 
that the Coulomb potential seen by the n electrons is that of a point charge 
at the origin; since there is at least one of the outer balls that is a distance 
2L from the origin, and since that ball contains at least one electron, we 
can safely say that the Coulomb interaction of the outer electrons with the 
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nuclei, i.e., the sum of the term (3.14) and the last term in (3.23), is less 
than — [Z — N')/3L. (The reason we wrote Z — N' instead of n(Z — TV 7 ) is 
that we do not know the positions of the other n — 1 electrons; they could 
be very far away.) 

To summarize the situation thus far, we have a negative Coulomb attrac- 
tion of order CL^ 1 , where C is a fixed constant. We have various localiza- 
tion errors of order L -27 , R 2 and also (R /r')(L-2R y y (l + \log(AR )\). 
These latter terms can be made arbitrarily small compared to CL _1 if we 
choose 1 > 7 > 3/4 and 

L 27 " 1 » Ro » L 1/2 . (3.24) 

Finally, there are the terms (3.11) and (3.18 - 3.22) which involve expectation 
values of linear operators a* in <J>(n) and in ^f(N'). These are dangerous 
looking terms; on the face of it they appear to possibly be of order L -1 , but 
we can make them all effectively vanish! 

To eliminate these terms we can make an anti-unitary transformation on 
ty(N') (or else on <&(n), but not on both) that will not alter the energy of 
each subunit or alter the Coulomb interaction. This anti-unitary is simply to 
replace a* by —a* and, simultaneously use complex conjugation to change 
ty(N') to its complex conjugate ^f(N'). In addition we apply the unitary 
operator W = YliLi > where a 2 is the second Pauli matrix in the usual 
basis in which a 2 has purely imaginary elements and a 1 and a 3 are real. 

The effect of applying this anti-unitary is to replace (3.11) and (3.18 - 
3.22) by their negatives, whereas all other energy terms remain unchanged. 
Note that the anti-unitary when applied to ^f(N') changes the sign of one of 
the factors in (3.18 - 3.22) only. It changes the sign of (^f(N'),pi^f(N')) be- 
cause of the complex conjugation and it changes the sign of (^(N'),ai^(N')) 
because of complex conjugation and W. The terms (^f(N'), A(xi)^(N')) and 
(^(N'), B(xi)ty(N')) change their sign because of the change of sign of the 
a*'s. Thus, each of this terms can be negated, and one choice or the other 
will make the sum (but perhaps not each individual term) of (3.11) and (3.18 
- 3.22) non-positive. □ 



4 Localization Estimates for 'Free' Electrons and 
Photons 

The main result of this section is Theorem 4.3 which shows how to construct 
a state in which the 'free' electrons and the field are localized. This was used 



686 



EXISTENCE OF ATOMS IN NON-RELATIVISTIC QED 



in Part 1 of the proof of Theorem 3.1. Part 2 of the proof of Theorem 3.1 
also uses the part of Theorem 4.3 relating to the field localization. 

The proofs in this section rely, in part, on the commutator estimates of 
Sects. 5 and 6. The Hamiltonian for the n free electrons is given in (2.3). 

4.1 Localization of the Electrons, but not the Photons 

LEMMA 4.1 (localization of electrons). Fix a radius Ro > 0. Then 
there exist (not necessarily disjoint) balls B\,...B n in M 3 , each of radius 
Ro, and a normalized vector VP in the physical Hilbert space TL(n) such that 
the electronic part of \£ is supported in Sl{B\, . . . , B n ) and with an energy 

(ty,H°(n)^) < E°(n) + bn 2 R Q 2 , (4.1) 

where b = 2ir 2 is twice the lowest eigenvalue of the Dirichlet Laplacian in a 
ball of radius 1. 

Conjecture: The proof does not tell us the location of the n balls. If they 
happen to be distinct then we can replace n 2 by n in (4.1). We conjecture 
that the theorem can be generally improved in this way, i.e., n 2 — ► n, with, 
perhaps, a different value for b. 

PROOF: Let e = bn 2 R^ 2 /2 and let $ be a normalized approximate 
ground state with error at most e/2, i.e., and ($, H°$) < E°(n)+e/2. 

Let B denote the ball of radius Rq centered at the origin in M 3 and let \ be 
a normalized, nonnegative, infinitely differentiable function with support in 
B. Define the function G of X = (x\,...,x n ) and Y = (yi,...,y n ) by 

n 

G(X,Y)= l[x(xi-V*i) , (4-2) 

where S n is the symmetric group. Clearly G is a symmetric function of the 
X variables and of the Y variables and, therefore, G(X, Y)$ is a valid vector 
in the physical Hilbert space for each choice of Y. 

It is obvious that 

P(X) := [ G(X,Y) 2 dY (4.3) 

is simply n! times the permanent of the nxn hermitian, positive semidefinite 
matrix Mjj := j R3 x( x i ~ y)x( x j ~ y)dy. It is a general fact that such a 
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permanent is not less than the product of its diagonal elements; in fact it is 
not less than the product of the permanent of any principal (n — m) x (n — m) 
submatrix and the permanent of its m x m complement [11], so P(X) > n\. 
In particular, a fact that we shall use later is that for each i, P(X) > 
nlMiiMii = n\Mii, where Mj j is the permanent cofactor of Mj j, i.e., it 
is the permanent of the matrix in which the i th row and column is deleted 
from M. (In our case, the assertion is obvious since every matrix element 
M id > 0.) 

We define W(X,Y) = G(X, Y)P{X)~ 1 / 2 , and using this we define 

tfy ■= W(X,Y)$ . (4.4) 

Our \l/ will be be Vl/y (up to normalization) for a choice of Y to be determined 
shortly. Since P(X) > n\ the multiplier W(X,Y) is C~. 

We proceed analogously to Theorem 3.1 of [7]. Consider 

£(Y) := (y Y ,H (n)* Y ) - [E°(n) + e + bn 2 R 2 } (tf y , tfy) . (4.5) 

Our goal is to show that J £(Y)dY < for a suitable choice of %. This will 
prove that there is a set of Y's of positive measure such that ^y / and 
also (* y ,i7°(n)*y)/(^y, *y) < [E°(n) + e + bn 2 R 2 }, which is what we 
wish to prove. 

It is obvious, from (4.3) that / W(X, Y) 2 dY = 1 and so 

J (tfy, #y) = ($, $) = 1. (4.6) 
In a similar fashion one sees that 
J (*y, [a l^-x^+tf/^y^F = ($, [aj^ |x i -x i |" 1 +^/]^) • (4.7) 

i<j i<j 



Next, we compute 

||(V Xj + iA(xi))*Y\\ 2 = \\(V Xi W)n 2 + 5tf (<&, (V Xi W 2 ) • (V Xj + iA(xi))*) 

+ ||^(V iEi +iA(x i ))$|| 2 . (4.8) 

The middle term vanishes when we integrate over Y since J W(X, Y) 2 dY = 1 
and hence J V Xi W(X,Y) 2 dY = 0. The last term gives us the required 
contribution of the kinetic energy to J(^Y,H°( n )^ f Y)dY in (4.5), again 
using the fact that fW(X,Y) 2 dY = 1 . The first term is ($,Fj(X)$), 
where 

Fi(X) = J |v Xi {G(X,Y)P{X)- l l 2 ^ dY . (4.9) 
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Our proof is complete if we can show that Fi{X) < 3e/2n, which we shall 
do next. We start with 

v Xi [g{x,y)p(x)-v*) = p{x)-v*v Xi G{x,Y) 

- (l/2)G(X,Y)P(X)- 3 / 2 V Xi P(X) . (4.10) 

If we square this and integrate over Y we obtain (recalling that X7 Xi P(X) = 
2 / G(X, Y)V Xt G(X, Y)dY and / G(X, Y) 2 dY = P{X) ) 

FAX) = -^y J \V Xi G(X,Y)\ 2 dY - Jp L^\V Xi P(X)\ 2 . (4.11) 

We shall ignore the last term since it is negative. 
In order to compute V Xi G(X, Y) let us write 

n n 

G(X,Y) = J2x(xi-yj)»j(X',r) :=J2aj(X,Y), (4.12) 

3=1 3=1 

where 

H (X',Y')= I[x(xi-Vnl) , (4-13) 

and where SVi-i denotes the set of bijections of 1, . . , i, . . , n into 1, . . , j, . . ,n. 
Then, 

j \V x% G\ 2 dY = jrjr j V Xi a r V Xi a k dY < nf^ [ W^\ 2 dY 

J j=l k=l J 3=1 

= njrj \{V X ){xi-y j )\ 2 dy j j H (X' ,Y f ) 2 dY' = nC iti J \V X (x)\ 2 dx , 

3=1 

(4.14) 

where Y' = (y±, . . . , yj, . . . , y n ) and where C^j = j n{X', Y') 2 dY' equals 
(n — 1)1 times M^j, the cofactor of M^j in the permanent of M. However, 
P(X) > nlMi^Ma = nlMij [Ci,i/(n - 1)!], as explained before. There- 
fore, P -1 j \V Xi G\ 2 dY < n'j |Vx| 2 - The same inequality holds for any 
i = 1, . . . , n which gives us a factor of n 2 altogether. 

At this point we wish to choose \ to be the lowest Dirichlet eigenfunc- 
tion of —A in the ball B. That function is not in C£°(B), but it can be 
approximated by such a function so that the error in n 2 f | Vx| 2 is less than 
e/2. □ 
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4.2 Localization of Photons 

Our next task is to produce a state in which the photons are localized. A 
definition is needed first. We say that the electromagnetic field in a state 
$ is supported in a closed subset S C R 3 if each component a\ of the 
field satisfies 1 1 o-a (a?) 1 1 = for all x ^ S. To construct a localized state 
from any given state <3? we use the representation (2.15) of Fock space. We 
suppress the space and spin variables of the electrons for the moment. For 
a smooth cutoff function < j(y) < 1 define the localization operator J 
on Fock space in the following manner. is still given by (2.15) but the 
vector \ii,mi; ••• ;i n ,m n ) is changed to 

J\h,m i; ••• ;i n ,m n ) = —^= a *(jf ii r^--a*(jf in r"\0), (4.15) 

ym\\ ■ ■ ■ m n \ 

Clearly, J is a linear, self-adjoint operator and 

\\a X {y)Jn = (4.16) 

for all y that are outside the support of the function j(y). Note that J is a 
contraction, i.e., < ||3>|| for all 

To work effectively with the operator J the following commutation re- 
lations will be useful later. Since the electron variables are not relevant for 
the calculation, we suppress them here. 

LEMMA 4.2 (Commutation relations for J). For any f in the single 
photon Hilbert space L 2 (]R 3 ;C 2 ) we have (with a\(f ) = f a\{y)f{y)dy, as 
before) that 

a x (f)J = Jax(jf), Ja* x (f) = a* x (jf)J (4.17) 

[ax(f),J] = Ja x ((j - 1)/), [a*x{f),J] = -a* x ((j - l)f)J . (4.18) 



PROOF: For any state in the Fock space we have that 
Mf)J*] n (yi,\i;...;,y n \ n ) 

n ,. 

= Vn + 1 Y[ j(vk) / f(y)j(y) [*] n +i (yA;yiAr,- ■ -;yn,K)dy 

k=l J 

= [Ja\Uf)*] n (viAi;---;y n ,K) ■ (4.19) 

All the relations follow immediately from this. □ 
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4.3 Localization of the Photons and the Electrons Together 

One would like to think that most of the photons ought to be localized near 
the electrons. This will be the case provided one replaces the state ^ of 
Lemma 4.1 by the ground state for the Hamiltonian H°(n) restricted to the 
states that vanish outside the set f2. Moreover, this state will have an energy 
close to the energy E°(n). The following theorem makes this precise. 

THEOREM 4.3 (Localized photons and free electrons). Fix radii 
Ro > and L > 2Rq. Then there exist (not necessarily disjoint) balls 
Bi,...B n in R 3 , each of radius Rq and a normalized vector <E>(n) in the 
physical Hilbert space TL{n) such that the electronic part o/3>(n) is supported 
in £l(Bi, . . . , B n ) and the electromagnetic field is supported in S = U" =1 Pj 
where Pi is a ball concentric with B, L but with radius L. 

The energy of (n) satisfies 

(*(n), fl°(n)*(n)) < E°(n) + b^ + (1 + | log(Ai? )|), 

(4.20) 

for any 7 < 1 and where b = 2ir 2 is twice the lowest eigenvalue of the 
Dirichlet Laplacian in a ball of radius 1. The constant c depends only on 7 
and is independent of Rq and L. 

PROOF: We start with the wave function given by Lemma 4.1. This 
fixes the balls B±, . . . B n , and hence the symmetrized product f2(-Bi, . . . B n ). 
The next step is to redefine the Hamiltonian H (n) by restricting the Hilbert 
space to the balls, i.e., we replace the space AL 2 (M 3 ;C 2 ) by the subspace 
of L 2 (Q; <g>™C 2 ) consisting of functions that are antisymmetric under the 
exchange of particle labels. (This makes sense because f2 is symmetric under 
exchange of particle coordinates.) The Laplacian is replaced by the Dirichlet 
Laplacian. 

A physical way to say this is that we add an infinite potential outside 
Q. This is not a sum of single particle potentials, but that is immaterial. 
The point is that by the methods of [7] there is a bound state, i.e., there is 
a state <3?£)(n) with lowest energy Ejj(n) (the letter D stands for 'Dirichlet') 
that satisfies Schrodinger's equation. (In fact, the methods of [7] are not 
needed to establish the existence of a ground state in this case since all finite 
energy states are evidently localized; this was noted earlier in [5], [1] and [8]. 
However, [7] is needed for the photon localization in the next step.) 

This ground state will obviously have a lower energy than the ^ given by 
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Lemma 4.1 since that $ automatically satisfies the Dirichlet boundary con- 
ditions, that is ($ D (n),H°(n)& D (n)) < (E°(n) + bn 2 R 2 ) ($ D (n), <f> D {n)). 

Next, we localize the photons in the set E. A standard IMS localization 
yields two smooth functions j\ (y) , ji (y) with 

ji(y) 2 +j 2 (y) 2 = i (4.21) 

with support of ji(y) in E, We also require that ji(y) is identically equal to 
1 on the set Uf =1 Qi where, for each i, Qi is a ball of radius L/2, concentric 
with Pi. Moreover we can assume that |Vjj(y)| < C/L for some constant 
C and i = 1, 2. We define J^u(n-) by using j\ in (4.15) and, with the help 
of (4.16), we use the localized state J'^rj(n), appropriately normalized, as 
a trial function. This function will be the required function $(n) of our 
theorem. 

The energy of the state J<&D{n) can be compared with the energy of 
$£>(n) by using the commutator formula 

(J* D {n), (H°(n) - E° D (n))J<Z> D (n)) = (j<S> D (n), [H°(n), J\* D {n)) . 

(4.22) 

An important point about having a ground state $o(n) is that one can 
derive infrared bounds for this state (see Sects. 5 and 6). All that is needed 
is that §d is a ground state, i.e., it satisfies the Schrddinger equation in 
order to apply the 'pull through formula'. 

Lemma 5.2 shows that the norm ||»7$d|| is close to one and Lemma 5.5 
shows that the right side of (4.22) is bounded as 

(j<f> D {n),[H°(n),J]<f> D {n)) 

~ (L- C 2R p (§) ( 1 + \^Ro)\)(J^D(n),J^ D (n)) , (4.23) 

for any 7 < 1, where the constant C depends on 7 but not on Rq and L. 
This shows that $(n) = J<S> D (n)/\\J<5> D (n)\\ satisfies (4.20). □ 



5 Commutator and Related Estimates 

In this section we prove various results stated in the previous sections, par- 
ticularly Lemma 5.5 which is used in the proof of Theorem 4.3. We also 
prove Lemma 5.6, which is not a commutator estimate; it is simpler. It is 
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needed to bound the terms (3.16) and (3.17), as we stated just after eq. 
(3.23). 

We will deal mainly with the Dirichlet ground state associated with elec- 
trons localized in the set Q(B\, . . . B n ). The bounds for the ground state 
describing electrons exponentially localized near the nuclei is easier and fol- 
lows in the same fashion. 

In this section and the next we denote the Dirichlet ground state 3>£>(n) 
simply by <E> in order to simplify the notation. 

Recall the definitions of in (4.21) and of the operator J in (4.15) 
which is defined by substituting j\ for j. An important operator in our 
analysis is the outer photon number, given by 

2 r 

Km = Yl / a l(y)a\(y)dy (5-i) 

x=1 Jj2(y)>o 
or, in terms of matrix elements, 

2 r 

(<&,AA out $) := / \\a x (y)n 2 dy (5.2) 
A=1 Jj 2 (y)>0 

We start with the following bound, which is a consequence of the infrared 
bounds proved in Sect. 6. It is used in the proofs of Lemmas 5.2, 5.3 and 
5.4. 

LEMMA 5.1 (Photon number is small far away from the electrons). 

For the Dirichlet ground state <J> of the free electrons localized in £l(Bi, . . . B n ) 
we have the bound ( with C independent of Rq and L ) 

($,AA$)<C(l + |log(A# )|) , (5.3) 
and for all 7 < 1, the bound 

(<l>,AA out <I>)<c(^ 2 ||<I>|| 2 (5.4) 

where the constant C depends on n, A, 7 but not on Rq an d L. Likewise, for 
the ground state ^ of the bound system given by the Hamiltonian H V (N') 
we have that ( with C independent of Rq and L ) 

(*,A/"#) < C , (5.5) 

and for all 7 < 1 and L > 2R$ that 

(*,AT ut*) < C (jL^J \\n 2 , (5.6) 
where the constant C depends only on 7. 
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Our goal is to prove inequality (4.23), which is Lemma 5.5 of this section. 
To do so, we shall need the following three lemmas, in which J is defined 
with j\. Recall that < ji(y) < 1 and j\(y) = 1 for y 6 UILi Qi ( see * ne 
proof of Theorem 4.3). 

LEMMA 5.2. For the normalized ground state $ ; we have for all < 7 < 1 
that 

l~\\Jn 2 <c{^) 2 (5.7) 

where C is a constant that depends on 7, n, A but not on Rq and L. Moreover, 
for an arbitrary state 

(J*,Kut JV) < (*,A/- out tf) . (5.8) 
PROOF: Formula (5.8) is immediate from 

oo n n 

n=l 1=1 1=1 



Next, note that 



l-i 



-n^w = Eiiii 2 (^ 2 fe) 



(5.10) 



k=l 



1=1 k=l 



(by definition the empty product equals 1). This is proved by inserting 
32 (Vn) = 1 — 3\{yn) on the left side of (5.10) and then repeating the process 
inductively. In particular, we have that 



U 11 



(5.11) 



k=i 



from which we obtain 



1 - \\J<$>f < ($,A/- out c&) < C {^j 2 (5.12) 



by Lemma 5.1. 



□ 



LEMMA 5.3. For the ground state $ and for every L > 2R$ we have that 



<c 1 



i=l 



(L-2i? )7 \L-yJ 



^Vl + |log(Ai?o)|)||J$|| 2 , (5.13) 
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for all 7 < 1. C is a constant that depends on 7, n, A but not on Rq and 
L. (Note that it makes no difference whether we use (pi + A(xi)) 2 or use 
its normal ordering since the commutator of ai and a* is proportional to the 
identity operator, which commutes with J '.) 



PROOF: We first calculate the commutator of (p + A(x)) 2 with J. 

n 

^r[( Pl +A( Xl )) 2 ,j] = 
i=i 

n 

(2pi ■ (K J] + [a*, J]) + [am, J] + [a* a*, J] + 2[a*a t ,J}) , (5.14) 



i=l 



where we abbreviated a\(h ] x (xi — •)) by and likewise for a*. (Note that 
the index j £ {1,2,3} of the coupling function is unimportant and will be 
suppressed from now on.) 

Step 1. The term Ya=i 2 Pi ' [di,J]$) is bounded, by Schwarz's 
inequality, by 

/ n \ l / 2 /n \ 1 / 2 

2lj2\\ P an 2 ) n2\\[ai,j]n 2 \ . (5.15) 

The first factor can be estimated simply in terms of the energy while the 
second factor will deliver the necessary decay in L. Using Lemma 4.2 the 
problem is reduced to estimating, for each fixed X = (x±, . . . ,x n ) G £1 and 
each i (with || • || denoting the norm in Fock space only) 

\\Ja(h(x i --)(j 1 (-)-l)M 

< \\a(h( Xi - -)Oi(-) - 1))*|| < /(! - h{y)) \\h(xi - y)a{y)ndy 

< f{l-h(y)) \h(x t -y)\\\a(y)<5>\\dy 

< (/ (i-ii(y)) |/^- y )|V) 1/2 (| (l-j!^))^^)^!! 2 ^ 172 . 

(5.16) 

The first factor in (5.16) can be bounded, using Lemma B.l, by 

1 \ 1/2 C 



{S {1 ~ JI{V)) (L-\R^ Xt ~ V?1HXt ~ v)?dV ) ~ 



(L - 2R p 

(5.17) 
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since, whenever 1 — ji(y) 7^ 0, the distance between X{ and y is at least 
d = (L/2) — Ro, by construction. If we now square (5.16) and integrate over 
X we get the desired decay estimate for (5.15). 

For the second factor in (5.16) we note that 1 — ji(y) < 32{y) 2 - This, 
together with Lemma 5.1, yields 

it wo* - •)(*(■) - d)*ii 2 ) ' < jrSm^ (§) < 5 - 18 > 

for all 7 < 1. 

Next, (J$,2E" = iK'K.^] $ ) = -2£™ =1 ([a;, Jjj^^andthiscan 
be estimated in the same fashion as before except that the estimate is in 
terms of 

instead of (<&, AT out &). On account of Lemma 5.2, this is bounded by 

($,A/- out $) . 

Hence, we obtain the same kind of bound as in (5.18), i.e., for all 7 < 1, 

lg(^, 2pi . K .,^)l< ir -|^(f) . 

Step 2. Returning to (5.14) we concentrate on the term [aiai^J] which 
can be written as ai[ai,J] + \ai,J]a,i. Using Schwarz's inequality 

{J^aiiai, J]$) = (a*J<f>,[ ai , J]$) < \\a*J<S>\\ ||K . (5.19) 

The second factor is treated in precisely the same fashion as in Step 1. The 
first factor cannot be estimated directly in terms of the energy, since the 
function is not an eigenfunction. This will be dealt with below where 
we estimate the term 1 1 a,* ^ 1 1 - 

The term 

(JQ,[ai,J]ai*) , (5.20) 
can be written, using Lemma 4.2, as 



(J$,Ja(h(xi - -)O'i(-) " l))a(h(x, t - ■))<&) 

= (j 2 ^,a(h( Xi - •))a(/i(x l - -)O'i(-) - 1))$) (5.21) 



696 



EXISTENCE OF ATOMS IN NON-RELATIVISTIC QED 



which, again using Schwarz's inequality, can be bounded by 

\\a*(h(xi - -))J 2 n \\a(h{ Xi - -)Oi(-) - . (5.22) 

As before, the first factor cannot be estimated in terms of the energy, since 
J 2 <& is not an eigenstate. Note, however, that 

\\ a *ji<f>f = \\ a J 2 <Z>\\ 2 + H/ifll^H 2 • (5.23) 

and the first term on the right side can be estimated by 

\\h\\ 2 (j 2 $,A/\7 2 $) < INI 2 ($,AA$) . (5.24) 

This follows from the formula 

oo 

^ = 5>*(/>(/i) > (5-25) 

3=1 

which is valid for any orthonormal basis {fj}, in which we pick fi(y) = 
h(x — y)/\\h(x — -)||, and from (5.9). (Here h is an abbreviation for the 
coupling functions.) Thus, using Lemma 5.1, 

| £ [wuJm I < C {L _\ Ro p (jj) ( X + I lo s( A ^o)l) • (5.26) 

Step 3. By taking adjoints the third term in (2), YH=i\ a *i a *i >«^] leads to 
the expression 

n 

-J2([aiai,J]J^) (5.27) 

i=i 

and can be dealt with in the same fashion as in Step 2. It remains to analyze 
Y!i=i[ a i a ii J] = YJi=i a *\ a ii J] + [a>*,J]<H- The fi rst term is estimated using 

(diJS, [<h, J]$) < \\aiJ$\\ || [a t , J]$|| , (5.28) 

while the second one can be written as 

-([ai,J\J*,a&) (5.29) 

which, once more by Schwarz's inequality, can be bounded by 

\\[(H, J]J$\\ \\ai*\\ . (5.30) 

Both these terms have been estimated previously. □ 
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We come now to the third lemma needed for the proof of Lemma 5.5. 
This lemma concerns only the real part of a commutator expectation value 
(5.31), but this is all we need for Lemma 5.5. The reason is that the total 
commutator in Lemma 5.5 is manifestly real, since <3?(n) is an eigenstate of 
H°(n) and J is selfadjoint. On the other hand, the piece of the commutator 
considered in Lemma 5.3 is also manifestly real and the only other part of 
H°(n) to be considered is the potential energy terms, which commute with 
J. Therefore, the commutator expectation value in (5.31) is, in fact, real. 
The proof of Lemma 5.4 is greatly simplified, however, by being able to 
ignore the (non-existent) imaginary part. 

LEMMA 5.4 (Commutator of J with the field energy). The ground 
state $ satisfies the bound 

R (Jd>, [H f , J]$)<j (1 + | log(Ai? )|) 1/2 (5.31) 

for any 7 < 1 . C is a constant that depends on 7, n, A but not on Ro and 
L. 

PROOF: It is convenient to write the field energy of a state 3> in the 
form 

($, H f $) = (2vr) 3 (oa(-)*, ^ « A (0^) = 

A 

^ J \x-y\ 

(see [10, Eq. 7.12(4)] ). Next, we note that the commutator expression 
(5.31) is given by 

First, we investigate the numerator of the sum of the two integrands, 
which is 

(a x (x)J$, a x {x)J<S>) - K (a x (x)J 2 $, a x (x)$) 

- H (a A (ar) J*, a A (y)J$) + H (a A (x)J 2 ct>, a A (y)<D) (5.34) 
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plus the same thing with y and x exchanged. For brevity's sake we have 
omitted the sums over A here and in the following. We note that the x-x 
terms (and likewise the y-y terms) cancel. This follows from 

(a x (x) J<&, a x (x) J<5) - 5ft {a x (x)J 2 <5>, a A (x)d>) = 
(a x (x)J$, a x {x)J<S>) - 5ft (Ja x (x)J<!>, a A (x)$) - 5ft ([a A (x), J] J$, a x (x)$) 
= »(aA(x)J*,[a A (x),J]*)-»([a A (x),J] J$,a A (x)$) , (5.35) 

which, together with Lemma 4.2, yields 

Cn(:r) - l)»(o A (x)J$, Ja A (a;)*) - (ji(x) - l)M{Ja x {x)J<f>,a x {x)<S>) 

= Cn(x) - 1) [St(ax(x)J$,Ja x (x)$) - ft(a x (x)J<S>,Ja x (x)<f>)\ = . 

(5.36) 



Now we deal with last two terms in (5.34). 

5ft (a A (x).7 2 <I>, a A (y)$) - 5ft {a x (x)J$, a x (y)J<5>) 

= 5ft ([a x (x), J] J$, a x (y)$) + 5ft (a x (x) JQ, Ja x (y)$) 
-$t(ax(x)J<S>,a x (y)J<5>) 
= 5ft {[a x (x),J] JQ, a x (y)$) + 5ft (a x {x)J$, [J, a x (y)} $) . (5.37) 

Again, by Lemma 4.2, this equals 

(h(x) - l)X(a x (x)J$,Ja x (y)*) - (h{y) - 1)» (a A (x).7$, Ja x {y)^) . 

(5.38) 

Commuting the J once more with the a's leads to 

tfi(x) - l)5ft ([ax(x), J] Ja A (y)<I>) - (j^y) - l)5ft ([a A (x), J] Ja A (y)$) 

+ {iiW-ii(y)}»(J r a A (x)$, l 7a A (y)$) . (5.39) 

We do not have to worry about the last term, since it is the product of a 
symmetric and an antisymmetric term in the variables x and y, and thus its 
x-y integral with \x — y|~ 4 in (7) vanishes. The other term, using Lemma 
4.2, is of the form 

{(ji(x) - If - (h(y) - l)(ji(x) - l)}5ft(Ja A (x)$,Ja A (y)$) . (5.40) 

Taking into account the terms with x and y exchanged we find that (5.33) 
equals 

4. y, J tow-M)?*^*)*.^)*)^ (M1) 
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Now we use Lemma 5.1 to get an estimate on the size of (5.41). Recall 
that the function ji(x), which defines J and which is defined in (4.21), is 
identically equal to 1 on U" =1 Qj, where Qi is the ball of radius L/2 equi- 
centered with the ball Bi. Moreover ji(x) = whenever the distance of x 
to the center of every Bi exceeds L. Write 1 = xi + X2 + X3 where xi 1S 
the characteristic function of U" =1 Qi and X2 is the characteristic function of 
the shell between U™ =1 Qi and £ = Ui =1 Pj. Finally, X3 1S the characteristic 
function of the outside region (on which j2(x) = 1 and ji{x) = 0). We note, 
for later use, that J %i < CnL 3 and J X2 < CnL 3 , where C is a universal 
constant. 



Next, we analyze each of the terms 



(jiW-hiyrfx^XjiyWi Jaxix)®, Ja x {y)^) 
' "' ' \x-y\' 



(5.42) 

Clearly, T± t i = T^3 = 0. To bound the other Tjj's we recall that 

(ji(x)-ji(y)) 2 <^\x-y\ 2 . (5.43) 
With this and Schwarz's inequality, Tjj is bounded by 

Cf MznJ^m x,(y)\\J a >(<M\\ 
j_j j \x y\ 

L 2 J \x-y\ l 

Denote ||aA(x)$|| by fix). Consider the terms i = 3 and j = 1,2. Using 
the Hardy-Littlewood-Sobolev inequality (see [10, Theorem 4.3] with 1/2 + 
5/6 + 2/3 = 2), we get the bound 

r 3J <-^||x3/||2||x,-/ll6/5, .7 = 1,2. (5.45) 
By Holder's inequality 

WXjfk/^WxMfh^Cn^LWfh, (5.46) 
and hence, for j = 1,2, 

T 3J < ^j-i®,^) 1 / 2 ^,^®) 1 / 2 . (5.47) 

Note that HX3/H2 is proportional to ($, A/"o U t^) 1//2 while ||/||2 is proportional 
to ($,AT$) 1 /2. 
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The next term to consider is i = 2 and j = 1, 2, 3. Again, H-L-S leads to 
the bound 

T 2J < ^11x^/11211x2/116/5 , j = 1,2,3. (5.48) 
Applying Holder's inequality yields 

||X2/||6/5 < IIX2||3 IIX2/II2 < Cn 1 / 3 L\\ X 2fh , (5-49) 
and hence the term with i = 2 and j = 1, 2, 3 is also bounded above by 

t ^-^^) 1 ' 2 ^,N ^) 1/2 , (5.50) 

where we used ($, A/" ut^ > ) 1 ^ 2 > ||X2/||2- It is the term (<3?,AA<1?) which yields 
the logarithmic term in formula (5.31). 

Finally, the term i = 1, j = 2 is the same as i = 2, j = 1, and the term 
i = 1, j = 3 is the same as i = 3, j = 1, both of which have been already 
treated. 

Collecting the estimates we have shown that 

Cn 1 / 3 

dt(J<Z>,[H f ,J} d») < i^($,AT$) 1 /2($ )A r out $)V2 

which, by Lemma 5.1, proves the lemma. □ 

We are now ready to prove the estimate stated at the end of the proof 
of Theorem 4.3, i.e., 

LEMMA 5.5. For all L > 2R we have the estimate 

(WM^,^)) <C/«A (1 + |]og(Afl0)|) 



(J$(n),J$(n)) - (L-2i? ) 7 V L7 
/or a// 7 < 1. T/ie constant C depends on n, A, 7 not on Rq and L. 

PROOF: The lemma is a direct consequence of Lemmas 5.2, 5.3 and 5.4. 

□ 

LEMMA 5.6 (Bound on the error terms). For each fixed x G IR 3 , 

for any 7 < 1, where D is the distance of x to the support of j±. The same 
estimate holds for a*(h{x — -)) 2 and for a*(h{x — -))a(h(x — •)) in place of 
a{h{x--)) 2 . 
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PROOF: Using Lemma 4.2, 

(J*, a(/i(x - -)) 2 ^) = (J$, Ja(h(x - -)ii(-)) 2 ^) = 

(a*(h(x - -)ii(-))^ 2 *,a(^ - -)j'i(-))*)- (5-53) 
By Schwarz's inequality this is bounded above by 

\\(a*(h(x - -)h{-))J 2 n Hh(x - • (5-54) 

Similar to the proof of Lemma 5.3, the second factor in (5.54) can be bounded 
as follows. 

\\a(h(x - < j h(y) \\h(x - y)a(y)*\\dy 

< J ji(y) \h(x-y)\\\a(y)*\\dy 

< (| ji(yf \h(x-y)\ 2 dy^j ' (| ||a(y)*|| 2 dy) 



(5.55) 

The second factor here is ($, N<&) 1 / 2 while the first factor can be estimated, 
using the fact that the support of j\ and the point x are a distance D apart, 
as 



^ (| n(yf \x - y\*\h{x - y)\ 2 dy^ < ^ (5.56) 

by Lemma B.l. 

The first factor in (5.54) can be bounded as follows. 

\\(a*(h(x - -)ii(-))^ 2 ^H 2 = (J 2 *,a(h(x - -)j(-))a*(h(x - -)ji(-))J 2 ^ 

(5.57) 

= {JH, [a(h(x - ■)j 1 (-)),a*(h(x - ■)ji(-))]J 2 ^ + \Hh(x - -) 3l {-))J 2 n 2 . 

(5.58) 

The first (commutator) term in (5.58) equals j \h(x — y)\ 2 ji(y) 2 dy || l 7 2( J > || 2 , 
which is bounded by C/D 2j . The second term was treated in (5.55) and is 
thus bounded by C/D 2 ~t($>,N$>). (Note that ($,A/"$) > (JQ,AfJ&).) 

It is immediate that similar estimates hold when aa is replaced by a* a* 
or by a*a. □ 



6 Infrared Bounds 



In this section we prove Lemma 5.1. It will follow from infrared bounds 
similar to the ones proved in [1] and [7], which have been proved to hold for 
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the electrons bound to the Coulomb potential. Those infrared bounds are 
not sufficient, however, for the localized wave function of the 'free' electrons. 
The chief reason for this insufficiency is that we need to know the dependence 
of the constants in the infrared bounds on the parameter Rq. Ultimately, 
the trouble stems from the fact that we do not know the positions of the 
n localized electrons, not even remotely. Thus, a direct application of the 
estimates in [7] would lead to constants that can grow conceivably as Rq or 
even faster for large Rq. This problem does not occur for the the bound 
electrons since, in that case, the electrons are localized by the Coulomb 
potential and the infrared bounds do not depend on the parameter Rq. The 
theorem below holds for the localized electrons, as well as for the bound 
electrons. The proof for bound electrons is easier and is omitted. 

As in Sect. 5 the Dirichlet ground state <&d(tz) for the 'free' electrons 
localized in in Q(B\, . . . , B n ) is denoted simply by <E>. Its energy is E^(n). 
The ground state for the bound system with Hamiltonian H V (N') is denoted 
by 

LEMMA 6.1 (Infrared bounds). The following infrared bounds hold for 

WMvn < i^xa^) , (6-i) 

||Sa(*)*|| < ^XA(k) • (6.2) 

The vector a\(k)& is a sum of n terms of the form e~ ik ' Yj Tj^\{k) where 
Tj t \(k) is given by (6.25) and satisfies the estimates 

The vectors Yj are defined below. The constant C depends on n, on the 
ultraviolet cutoff A and it is a monotone decreasing function of Rq . Similar 
bounds but without the factor Rq hold for the bound electron ground state ^ . 

PROOF: We prove first the bound (6.2) in detail. As a first step one per- 
forms an operator valued gauge transformation (see [7, Eq. 47] ) by applying 
the unitary operator (in which A is the vector potential (2.4)) 

m 

U{x) = eoq>[-iVa J^-(:r)(a; - Yj) ■ A(Yj)] (6.4) 

3=1 

to the wave function in each of its variables, i.e., 

n 
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Here, (f>j is a suitably chosen smooth function of compact support and the Yj 
are suitably chosen vectors. Note that the factors in the product commute 
since A(x) commutes with A(y) for all x and y. 

Next, we describe the functions (f>j. Consider the balls Bi(2Ro), which 
are concentric with the Bi but with twice the radius, and group them into 
clusters according to whether they overlap or not. We denote the number of 
these clusters by m. Such a cluster of balls Cj has a diameter that is bounded 
above by 4Rq times the number of balls in the cluster and hence bounded 
by AuRq. Denote by Yj the center of the cluster Cj which is defined to be 
the center of the smallest ball that contains all the balls Bi(2Ro) belonging 
to that cluster. We choose functions that are smooth , supported in 
the union of the balls Bi{2Ro) that belong to the cluster Cj and that are 
identically one on the union of the balls Bi. In particular ^j=i<A?( x ) = 1 
for x e Uf =1 Bi. 

The gauge transformation IA transforms the Hamiltonian H°(n) into the 
Hamiltonian 



H°(n) = UH°(n)U* = £ [( Pi + ^A{x t )) 2 + 9 -^a t - B( Xi j 

i=i 

+ a Y.T^-\+Hf (6-6) 

i<j \ Xl x i\ 

where the new field A(x) = UA{x)U* + a'^UpU* = A(x) + a'^UpU* is 
given by 

m m 

A(x) = A(x) -J2<t>j(x)MYj) - V^-(x)(x - Yj) ■ A(Y 3 ) , (6.7) 

3=1 3=1 

in which A(x) is still given by (2.4). The transformed field energy Hf is 
given as in (2.9) but with a\(k) replaced by the transformed creation and 
destruction operators 

b x (k,X) =Ua x (k)U* = a x {k) - i^aw x {k,X) (6.8) 

with 

Wx (k,X)=Y / EM^-y3)- £ T^e-^Xx(k) ■ (6.9) 
i=i j=i 

As before, the letter X denotes the vector (xi, . . . , x n ). We note that 

a x (k)<5> =U* \a x {k)$ -i^w x (k,X)$] . (6.10) 



704 



EXISTENCE OF ATOMS IN NON-RELATIVISTIC QED 



Since <& satisfies the Schrodinger equation we can apply the standard 
pull-through formula [1, 7] and compute 



H°(n)-E° D (n) a A (fc)$ = H°(n), a A (fc) U = 



2xA(k)V^\k\' 1/2 ex(k)-^2( Pl + ^A(x l )) I J>~*^ - e^^fo) \ d> 

i=l I j=l 



n m 



— y/\k\ 



ex(k) 

i=l j=l 

+ 4xA(k)^ kAe ^ .J2<7je- ik Xi *- \k\bx(k,X)* . (6.11) 

2 v\ k \ 



The term Yl^Li e lk Xi( t>j( x i) stems from the commutator of a\(k) with 
A(xi), which yields a term proportional to e~ lk ' Xi without any summation 
on j. Using, however, the relation J2Y=i <Aj'( x *) = which is valid for any 
Xi in Uj =1 Bj, we get that e~ lk ' Xi = Y^JLi e ~ %k ' Xi 4>j{ x i)- Note that X{ has to 
be in U™ =1 Bj for, otherwise, $ vanishes. 

Likewise, the second term in the above formula appears to have bad 
infrared behavior, but in order that does not vanish it is necessary that 
be in one of the balls. But then V(j)j(xi) = since <pj is constant, and this 
term does not contribute. 

Since H°(n) — E Q D (n) is nonnegative the operator H°(n) — E® D (n) + \k\ 
has a bounded inverse R(k). Thus (6.11) leads to the equation 



S A (fc)cD = ^£e-^2xAW |fcr 1/2 £A(fc) • R(k) 

n 

x J> + yfaA( Xi ))(l - e lk < Y ^)Mx^ 
i=i 

, 1 



+ ^^^XA(fc) fcA ^ fc) • R{k)aie~ ik ' Xi § + ^|fc|#(£>A(fc, X)* . 

(6.12) 
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As in [7], simple estimates, using Schwarz's inequality, lead to 

n 

\a x (k)n < 2xA(k)V^ \k\ 1/2 \\R(k) + ^A{ Xi )) 2 R{k)\\ 1/2 

i=l 

£(X>i-^(*i)*r 2 ' 

i=l 3=1 

+ ^W)v^ n\k\V 2 \\R{kj*\\ + \k\\\R{k)\\\\w x {k,X)l\\ . (6.13) 

Note that the index j in the first summand is determined by the ball to 
which the electron i belongs. Therefore, \xi — Yj\< 2uRq. 

Lemma A.l states that 

i— n 

9 -^Y,^-B(x J ) + H f + C>0 (6.14) 

3=1 

where 

c = ^2^ 2 aN 2 J XA(k) 2 dk . (6.15) 

t 



Thus we also have that 

N 

Oj ■ B( Xj ) + H f + C > , (6.16) 



2 

j'=i 



(note that the gauge transformation U commutes with B{x)) and hence 

n 

(5 : = \\R(k) Y J (P l + V^A(xi)) 2 R(k)\\ 
i=i 

< \\R(k) ^[(Pi + V^Mx t )f + ^a t B( Xi )} +H f + C ) j R(k)\\ 

= \\R{k)(H°(n) + C^R{k)\\ . (6.17) 
Thus, by subtracting and adding E® D — \k\, 

< \\R(k) (H\n) - E%{n) + \k\) R(k)\\ + (\E° D (n) - \k\\ + C)\\R(k) 2 \\ 
= \\R(k)\\ + (\E° D (n)-\k\\ + C)\\R(k) 2 \\ (6.18) 

where the constant C depends on A and n. For the last term in (6.13) we 
have that 

m n 

\\w x (k,X)§\\ <Mk)\k\~ 1/2 t252\\\ x i- Y j\M*i)*\\ ■ ( 6 - 19 ) 

j=i i=i 
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Since \\R(k)\\ = l/\k\, we have that ||fl(fc)$|| < l/|ife|. 

By combining these estimates with (6.13) and using \xi — Yj\(f)j(xi) < 
2nRo, we obtain the bound 

\\ax(k)n <C J R |A:r 1/2 ||$||xA(fc) , (6.20) 

where the constant C depends on A, n and the energy. This estimate carries 
over to the state by (6.10) since w\(k,X) applied to $ satisfies the same 
estimate, as we see in (6.19). Note that the energy Ep(n) does depend on Rq 
but it is monotone decreasing as a function of i?o (by the variational principle 
for Dirichlet boundary conditions) and it is uniformly bounded below. 

Next, we observe in (6.12) that a\(k)& is a sum of m terms of the form 
e~ tk ' Y JSj(k) where 

S jtX (k) = ^XA(k) \k\- l ' 2 e x (k) ■ R(k) 

n 



n 

X 

i=l 



2 V\ k \ T^i 

n 

+ i^XK{k)\k\ l ' 2 ex{k) ■ R(k) M x i)( x i ~ Y J)* > ( 6 - 21 ) 

i=i 

where we have used the identity 

n m n 

J2^e~ ik ' x ^ = j2 e ~ ik ' Yj J2 ai( t ) ^ x ^ e ~ ikiXi ~ Yj) ^ ■ ( 6 - 22 ) 

1=1 j=l 1=1 

Since by (6.8) 

a x (k)$ = W [a x (k) - i^w x {k, X)] $ , (6.23) 

we obtain that 

n 

a x (k)* = J2 e ~ ik ' Yj Tj,x(k) (6-24) 

3=1 

where 

f j7X (k) = W % x (k)$ - i^w x (k, X)} $ . (6.25) 



Differentiating these expressions with respect to k and proceeding in the 
same fashion as in the proof of (6.2) yields the estimate (6.3). 
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Differentiating the polarization vectors (2.7) yields the factor \Jk\ + k\ 
in the denominator of (6.3) . The details of the calculation are the same as 
the ones in [7] and are omitted. The bound (6.1 is considerable easier, since 
its proof does not require the gauge transformation. Otherwise the proof is 
word for word as the one above. Finally, the proof of the infrared bounds 
for ^ is a word for word translation of the one given in [7]. Note that the 
localization radius does not show up in this calculation since the electrons 
are exponentially localized in the vicinity of the origin. □ 

Finally, we come to the main application of the infrared bounds proved 
in this section. 

PROOF OF LEMMA 5.1: Using (6.1) we can write 
($,M,ut$) = / \\a{x)n 2 dx= [ \\y2T jtX (x-Yj)\\ 2 dx (6.26) 

Jj2(x)>0 Jj2(x)>0 j=1 

which, by Schwarz's inequality is bounded above by 

WT^x-Y^fdx . (6.27) 

r)>0 

For x in the support of ji we have that \x — Yj\ > L and hence 

/ \\Tj,x(x - Y^fdx < -L / \x\^\\T hX (x)\\ 2 dx . (6.28) 

Jj 2 (x)>0 L^ J 

This last term can be related to the derivative in k space of the function 
Tj t x(k) by the formula 

f „ „ r(v k T h x{k)^ k T h x{k')) 
j \x\^\\T hX (x)\\ 2 dx = C,j ± \ k _ k i\2 1+ i ~ dkdk ( 6 - 29 ) 

where the constant C 7 is given by 

C, = (^r 3/2 ^^y (6-30) 
Indeed, writing Qj \{x) = xTj iX (x) this formula follows from 

y(k),Q jA (k')) 



n „ 



J \x\*r- 2 \\Q j>x (x)\\ 2 dx = C, J V | fc _ fc/[27+1 J dkdk' (6.31) 
[10, Corollary 5.10] and the fact that 

(QjAk^QjAk')) = (v k T j:X (k),V k T j:X (k')) . (6.32) 
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Using the bound (6.3) a straightforward calculation shows that the func- 
tion ||VfcT} ) A(fe)|| is in L p for all p < 2. (The relevant term in (6.3) is the 
first term on the right side.) Using Schwarz's inequality and the Hardy- 
Littlewood-Sobolev inequality [10, Theorem 4.3] we can therefore bound 
(6.29) by 

2/p 



\V k T hX {k)fdk 



< CR , 



(6.33) 



with p = 6/(5 — 27), which is strictly less than 2 for 7 < 1. To prove (5.3) 
we write for some < H 

($,AA$) = V/ \\a x (k)$\\ 2 dk + Y] [ \\a x (k)<S>\\ 2 dk , (6.34) 

x J\k\<H ^JH<\k\ 



\k\ 

and, using (6.2) and (6.1) we get 

CRo [ 777 X/vdk + C [ T^rXhdk , 

J\k\<H \ k \ JH<\k\ \ k r 

which, optimized over H, leads to (5.3). 



(6.35) 



The proof for the state * is carried out in precisely the same fashion. □ 



Appendix A 



LEMMA A.l. On A N L 2 (R 3 ; C 2 ) <g> T we have that 



N 

£ °j ■ B{ Xj ) +H f + ^9 2 aN 2 I XK{k) 2 dk > . (A.l) 
.7=1 



PROOF: The magnetic field operator can be written in the form 

2 r 

£ / [c*x(k)a x (k) + c x (k)a* x (k)]dk (A.2) 



A=l 



where 



With this notation we can write the Hamiltonian (A.l) as 

2 



(A.3) 



a x (k)®I+jj^c* x (k) 



E/i*i 

\=i j 

2 r 1 
E J ^Wcx(k)dk 



a x (k) ® I + -^c x (k) 



dk 



(A.4) 
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Here, I denotes the identity operator on spin space. The first term is non- 
negative and crude estimates on the second yield the lemma. □ 



Appendix B 

In this section we prove the estimates on the coupling functions h\(y) which 
are defined by 

h\(y) = ±-J -J= £ \(k) XA (ky k -*dk . (B.l) 

It is important to choose the polarization vectors carefully in order that their 
Fourier transforms (from fc-space to y-space) have nice decay properties as 
\y\ tends to infinity. We shall express these decay properties in an integrated 
form. The reason for that is that the decay is not uniform with respect to 
the direction of the y variable. Recall the definitions (2.7). 

LEMMA B.l (Decay of the coupling functions). For any 7 < 1 there 
is a finite constant C (7) such that 

EE / \y\^\h\(y)\ 2 dy<C( 7 ) . (B.2) 

i=l A=l J 

PROOF: First we compute the gradient of h\ in fc-space. It is elementary 
that 

\V-^=e x (k)\ < 7= £3=5 , (B.3) 

where C is some constant. Because xa is smooth, 
3 2 

EE/ Iv^n^WxAWr^ < c( P ) (b.4) 

i=l\=l J vl fe l 
for any p < 2. We proceed as in (6.29), (6.30) and write 

EE/ \v\^mv)\ 2 dy = Cr EE / V }k-'jSli k>) dkdk ' • (B-5) 

i=l A=l i=l A=l 

Again, by the Hardy-Littlewood-Sobolev inequality [10, Theorem 4.3] this is 
bounded by 

2/p 



3 2 



EE/ \V^=e x (k) XA (kWdk 



=1 A=r 



< C p C(p) 2 /P , (B.6) 



where p = 6/(5 - 27) < 2 if 7 < 1. □ 
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